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Abstract A general framework that incorporates the Iroshnikov-Kraichnan (IK) 
and Goldreich-Sridhar (GS) phenomenalogies of magnetohydrodynamic (MHD) tur- 
bulence is developed. This affords a clarification of the regime of validity of the 
IK model and hence help resolve some controversies in this problem. This general 
formulation appears to have a certain robustness with respect to the inclusion of 
compressible effects. 



Magnetohydrodynamic (MHD) flows that occur naturally (like astrophysical sit- 
uations) and in modern technological systems (like fusion reactors) show turbulence. 
Early theoretical investigations of MHD turbulence considered the isotropic case. On 
the latter premise, Iroshnikov p] and Kraichnan [2] (IK) made arguments a la Kol- 
mogorov [3] and proposed that statistical properties of the small-scale components 
of the velocity and magnetic fields 

*are controlled by the shear Alfven wave dynamics; 

*show, in the limit of large viscous and magnetic Reynolds numbers, some 
universality in the inertial range; 

and gave for the total energy spectral density E[k), the behavior E(k) ~ k~ 2. Mont- 
gomery et al. [I] and [5], Goldreich and Sridhar [6] - [8] (GS) pointed out that the 
isotropy assumption in the IK theory is not a very sound one in the MHD case, 
thanks to the magnetic field of the large-scale eddies, and the GS theory [8] gave 

for the energy spectrum in the plane transverse to the magnetic field the behavior 

_ 5 

E(k±) ~ k ± 3 . However, EDQNM closure calculations of two-dimensional (2D) MHD 

turbulence dominated by non-local interactions in the presence of a mean magnetic 

field (Pouquet et al. [9]) and DNS of MHD turbulence in an applied magnetic field 

(Maron and Goldreich [10], Muller et al. [1 lj and [12]) showed that the transverse 

_ 3 

energy spectrum is close to the IK theory E(k±) ~ kj_ 2 . On the other hand, the 
3D DNS (Muller et al. [13], Haugen et al. [H]) of MHD turbulence, numerical cal- 
culations of 2D MHD turbulence dominated by local interactions (Fyfe et al. |15j) 
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and solar wind measurements (Leamon et al. [16J, Goldstein et al. [T7]) confirm the 
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GS spectrum E(k±) ~ k ± 3 . A resolution of this apparent conflict (Boldyrev [ [18J ) 
requires clarification of the regime of validity of the IK model which is the objective 
of this paper. A general framework that incorporates IK and GS phenomenalogies 
is developed to accomplish this objective. The IK and GS hypotheses can be shown 
to follow from the formal analogy between the hydrodynamic and MHD spectral 
energy density expressions. This general formulation is extended further to include 
compressibility effects. 

Let us write the spectral energy density E(k) as 

V 

E{k) 



k 2 T 

V being the characteristic velocity of the spectral element k. The hydrodynamic eddy 
turn-over time r given by 



then becomes 

1 



k?E2 



(3) 



(2) implies that the energy transfer in the hydrodynamic case is local in the spectral 
space which reflects the fact that a large-scale velocity field can be transformed away 
via Galilean invariance. 
If we use the relation 

T ~ — (4) 

e being the mean energy transfer rate, (3) leads to the Kolmogorov [19] spectrum 

E{k) ~ei*rf. (5) 
One may write for the MHD case, in analogy with (1) (Shivamoggi |20j), 

m ~ ^ ( 6) 

f being the MHD turn-over time, and Ca being the velocity of Alfven waves in the 
total magnetic field - 

C A = C Ao + C A (7) 
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where, 



Ca = ^: ) Ca = ^- (8) 

B is the applied magnetic field B is the magnetic field of the large-scale eddies, 
and k\\ and k± are the wave number components parallel and perpendicular to the 
magnetic field, k • Ca is a measure of the energy involved in the magnetic field-line 
bending (without which the magnetic field has no participation in the dynamics). 
Combining (6) with (3), we obtain^ for the MHD eddy turn-over time f - 



where, 



r~r(l + ^] (0) 
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k\\C Ao k±V k L C A . 



(9) implies that the energy transfer in the MHD case is non-local in the spectral space 
which reflects the fact that a large-scale magnetic field cannot be transformed away 
via Galilean invariance [2]. On the other hand, a mean magnetic field suppresses 
energy cascade parallel to it via more rapid Alfvenic decorrelation in the parallel 
direction so the energy spectrum would essentially be determined by 2D fluctuations 
[5]. We have from (9) 

T Z TA (Ha,b) 

(11a) corresponds to the case with a very strong applied magnetic field (B a 3> B) and 
represents the IK hypothesis, so the IK model is pertinent for a strongly magnetized 
plasma which is indeed anisotropic (invalidating the isotropy assumption in the IK 
model), (lib) corresponds to the case with a very weak applied magnetic field (B a 
B) and represents the GS hypothesis. 

(11a) leads to the IK spectrum - 



E(k ± ) ~ e^cik-J (12) 



-'^We assume that MHD turbulence is in a steady state so there is equi-partition of kinetic and 
magnetic energy associated with the advent of an Alfvenic state (V = ±B)(Matthaeus and Mont- 
gomery [2T]). 
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while (lib) leads to the GS spectrum - 
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E(k ± ) ~ £3k_J. (13) 

Let us now extend this formulation to include compressibility effects. It may 
be mentioned that the effects of compressibility on MHD turbulence are not well 
understood (Lithwick and Goldreich [24], Cho and Lazarian [25]). The following de- 
velopment should therefore be viewed as tentative. Assuming barotropic fluid and 
adiabatic flow processes, scale invariance of compressible Navier-Stokes equations 
leads to (Shivamoggi[22],[22]) ^ ne following scaling behavior of the velocity and den- 
sity increments - 

V(£) ~£ a ,p(£) (14) 

which implies the relation - 

2 

v\ ~ 



P-Po I- (15) 

c being the speed of sound and 7 the ratio of specific heats of the fluid. 
Writing the spectral energy density now as 

E(k) ~ ^ (16) 
the hydrodynamic eddy turn-over time r, given by (2) then becomes 



7-1 

2 7 



1— 7 1—37 



r ~ po~< c~^E^rk~. (17) 

In the zero-compressibility limit (7 — > 00), (17), of course, reduces to (3). Using 
the relation (4), (17) leads to the compressible energy spectrum ([22], [23]) - 

E(k) ~ pf^c'^e^k'^ . (18) 
One may now write for the MHD case, in analogy with (16), 

Combining (19) with (15) and (17), we obtain for the compressible MHD eddy 
turn-over time f c - 

f c ~ r M + y J (20) 
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which is identical to the incompressible MHD result (9)! Therefore, (9) is apparently 
form invariant with respect to the inclusion of compressibility effects. 

On the other hand, (20) implies that the IK spectrum in compressible MHD 
becomes 

E(k ± ) ~ pt^c-^e^C^k]^ 1 (21) 
while the GS spectrum in compressible MHD becomes 

E(k ± ) ~ pf^c'^e^k'^ . (22) 

(21) and (22) reduce to (12) and (13), respectively, in the zero compressibility 
limit (7 — > 00). 

It may be mentioned that Cho and Lazarian [25] found that the fast-mode turbu- 
lence, in the low-/? limit (which corresponds to the IK regime), exhibits the Zakharov- 
Sagdeev spectrum E{k) ~ k~^ (Zakharov and Sagdeev [26]). On the other hand, 
Lithwick and Goldreich [24] conjectured that, in the high-/? limit (which corresponds 
to the GS regime), the Alfvenic turbulence and slow-mode turbulence exhibit the 
scaling E(k) ~ k~ 3. 
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